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Down to the wire 


Part Two of a four-part series introducing the subject of graphics from a 
programming point of view. In this article John Knight, Paul Gordon and 
Philip Morris will develop the wireframe representation of our PCW object. 


W. will begin to produce some real 3D graphics 
in this part of the series. In the first article we 
introduced the subject of 3D graphics and the data 
structure we are going to use. Here we will discuss 
the matrix operations applied 
to the data structure to 
produce transforma- 
tions such as 
perspective, ro- 
tation, transla- 
f tion and 
“scaling. 
“et We will look 
at the subject of 
clipping to a viewing 
volume and the genera- 
tion of a wireframe image from our 
polygon/vertex database. 

The next part of this series will 
start to add the realism — the use of z- 
buffers and other hidden surface re- 
moval techniques will be explained. 
The Scan Converter method which 
can be enhanced to provide all the 3D 
effects we need to produce the fully 
rendered ‘PCW’ will be introduced 
next month. 


Matrix 
operations 
As we are using 
the computer as 
a visualisation 
tool we want to 
beabletomanipu- 
late any object de- 
scribed in its memory 

as we would a real ob- 
ject in our hands. The mo- 
tions traced out by ourhands 

can be described mathemati- 

cally as sets of linear transformations: i.e. transla- 
tion of the object in space, rotation about any axis, 
scaling and perspective. 


In general, for any single view of an object on the 
computer screen at least one of each of the 
aforementioned transformations must be applied 
to all points defining the polygons forming the 
object. A concise and economic way of applying 
the set of transformations is to use 2D matrices. 


Matrix representation 
A point in world space defined by its x, y and z 
coordinates can be written in a row matrix: 


ya 


An operation on this point can be written as a 
3x3 matrix, for example, scaling the point’s position 
by 2, thus: 


| pe eae 
Pee! ee 


When matrices are multiplied together, a single 
row matrix results containing the new point’s 
coordinates: 


ae 
(x y a(6 7 2 0) =ax By 229) 


In this way any rotation, scaling or shearing 
operation can be applied to a set of 3D points. 
However, in order to incorporate translation and 
perspective into the same framework we must use 
homogenous coordinates. 


Homogenous coordinates 
The homogenous coordinate matrix fora 3D position 
vector is written as 
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where h is the ‘distance’ in the fourth 
dimension, always set to 1 on input to 
the matrix. This coordinate can be oper- 
ated on using a 4x4 matrix, the upper 
left 3x3 component being identical to 
that described above, thus 


( Ketvis Yse Bewesl®) 


The multiplication of the homog- 
enous position vector and transforma- 
tion matrix produces the new coordinate 
in the same fashion as the 3D point (Fig 
1). 


SS 


ar 


As can be seen from the result, the terms in the 
bottom row of the transformation matrix allow 
translation in any of the coordinate axes. Each term 


H iiead « Sakae eed (ax + ey + iz + hm, 
e f eg hy. bx +fy+jz+hn, 
Bi Tokay i jf k 1 | cet+gy+kz+ho, 
Aleit dane shih dx + hy + Iz + hp) 
ran 
Nw 


in the transformation matrix, a through p, has a 
specific effect on the output coordinate. 
These are briefly described below and 
in more detail in the panel ‘Matrix Trans- 
forms’, page 284. 

The leading diagonal (a, f, k) { 


Fig 2 
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2 Rotate, 
matrix R1. 
3 Translate back to object origin, matrix 
T2; 

This can be achieved by multiplying all object 
coordinates by the single matrix V formed thus: 

v = T1 R1 12 

The order of the multiplication must be the order 
in which the transforms are applied to achieve the 
required result. 


Transform Matrix Matrix _multiply(mi1,m2) 


Scaling in x, y and z directions respec- 
tively. 

Bottom right (p) 

Global scaling as it effects only the h 
term; hence when projected to 3D space 
by dividing by h, the complete set of 
points is affected. 

Upper left 3x3 sub-matrix 

Facilitate application ofrotations about 


Transform Matrix R; 
int row,column; 


for(row=0; row<4; row++) 
for(column=0; column<4; column++) 
R.e[row] [column] = 0; 


for(row=0; row<4; row++) 
for(column=0; column<4; column++) 


any of the coordinate axes (see panel, { 

‘Matrix Transforms’). R.e[row] [column] += M1.e[row][0] * M2.e[0] [column]; 

Bottom row (m, n, 0) R.e[row] [column] += M1.e[row][1] * M2.e[1] [column]; 

Translation in x, y and z directions R.e[row] [column] += Ml.e[row] [2] * M2.e[2] [column]; 

respectively. R.e[row] [column] += Ml.e[row] [3] * M2.e[3] [column]; 
“™ Right column (d, h, J) } 


Non-zero values introduce perspective 
projections onto the x=0, y=O0 and z=0 +} 
planes respectively. Multiple instances 

of non-zero values produce two- and 
three-point perspective projections. 

The result of the transformation is a 4D homog- 
enous coordinate which must be projected back 
into 3D real space by dividing through by h. This 
projection allows the element of perspective to be 
introduced. 


Matrix concatenation 

Matrix representation provides us with asingle tool 
to apply any transformation to a 3D world space 
vertex. Once we have the set of transformations 
necessary to produce the required view of our scene 
on the computer screen in matrix form, we can use 
multiplication to produce a single matrix which 
applies all operations simultaneously! 

For example, to rotate an object about a given 
vertex in space v we must use a set of three transfor- 
mation matrices: 

1 Translate object so v is origin, matrix T1. 
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return R; 


This matrix combined transformation need be 
calculated only once per complete scene view, all 
points in the scene being operated on identically by 
the matrix. This method is memory efficient and 
fast. 

We can define two structures to represent a 
homogenous coordinate and the 4x4 homogenous 
transformation matrix: 


typedef struct 
{ 

float e[4] [4]; 
} Transform_Matrix; 


typedef struct 
{ 
float e[4]; 
} Homogenous_Coord; 
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Once set to contain the values for a single trans- 
formation, two structures oftype Transform_Matrix 
can be combined using the routine in Fig 2. 


can display on the screen, and by the end ofthis part 
we will be able to display a wireframe of our scene. 


Viewpoint 

Having defined the data structure that represents 
our scene, the next definition required is our view- 
point. For simplicity we consider the eye as a point 
light receptor, and this means we can define a point 
(in world coordinates) from which we are observ- 
ing the scene. 

When specifying a aaa we must not only 
consider the position of the eye, but also the ‘line of 
sight’: in other words, a line which joins the eye 
point to somewhere on the scene along which we 
are looking. 

The viewpoint and line of sight are normally 
defined together by a single ‘viewing vector’ which 
is a convenient mathematical (geometric) way of 
specifying a point and a direction. 

Having defined our viewing vector we must 
relate this to the screen coordinate system. The 
most convenient way of doing this is to assume that 
the viewing vector is along the screen z-axis, and 
that the eye is at the screen origin. This can be seen 
in Fig 3. This diagram shows a simplification often 
made for specifying the viewing vector: it is as- 
sumed that the line of sight is always from the eye 
point to the world coordinate origin. This means 
that the viewing vector can be defined simply by 


wy 


Fig 3 A suitable C function can be written to calculate | the eye point. 
the viewing transformation matrix for each input of 
scene viewing parameters, defined next. Specifying the Eye Point 
As mentioned before, the eye point is defined in 
3D Transformation world coordinates, and thus could be specified 
Last month we discussed how we could represent | using the same system as used for vertices in the 
a scene we wanted to render, and we talked about | scene. However, itis often desirable to view a scene 
vertices, polygons, beziers, coordinate systemsand | from different sides or above/below, and this posi- 
so on. tion can be more easily specified in polar coordi- 
Fig 4 This month we address the techniques of | nates (with rrepresenting the distance from the eye 
Ig 


converting our representation into something we 


to the world origin, T representing the angle from 


| Mat trans_eye(r, 
i float xr, 


theta, phi, 
theta, phi, 1_th, 


do thy 
1_ph; 


1_ph) 


/* generates a translation from world origin to eye origin */ 
{ 
float sin_y_rot, cos_y_rot, sin_x_rot, cos_x_rot; 
float r2; 


Vec eye=vec3d(); 
Mat transform = mat3d(); 


sin_y_rot = sin (-theta); 
cos_y_rot = cos (-theta); 
sin_x_rot = sin ((pi/2)-phi); 


cos_x_rot cos ((pi/2)-phi); 
S23 eo 
eye[1] 
eye [0] 
eye [2] 
eyelen=r; 


sin (phi); 
cos (phi); 
r2 * sin (theta); 
r2 * cos (theta); 
/* eye to origin distance */ 


transform=mtranslate3d (identity, -eye[0],-eye[1],-eye[2])i /* STEP.1 .*/ 
transform=mroty3d(transform, sin_y_rot, cos_y_rot); " /* STEP 2 */ 
transform=mrotx3d(transform, sin_x_rot, cos_x_rot); /*> STEP.:3 */ 


freevec3d (eye) ; 
return (transform) ; 
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the world z-axis and 9 repre- 
senting the angle from the 
world y-axis). In this way it is 
possible to view the scene from 
all sides simply by changing 
T, or to view it from above or 
below simply by changing 9, 
and so on. Fig 3 shows the 
polar eye point representation. 


The ‘Viewing 
Transform’ 

The first stage in the so-called 
‘viewing pipeline’ is the trans- 
formation of the scene into the 
screen coordinate system. We 
shall consider the case of a 
scene specified entirely with 
respect to a single world coor- 
dinate system — the use of 
object coordinate systems is 
simply an extension of this 
technique. 

To determine the transfor- 
mation required to map all 
vertices onto screen coordi- 
nates, it is sufficient to derive 
the transformation required to 


eee ccc 
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map the world coordinate origin onto the screen 
origin. This mapping can then be applied to all 
vertices in the world to produce a screen 
equivalent. 

If we assume the simplification 
described in the ‘Viewpoint’ sec- 
tion where we shall look along the 
screen Z-axis at the world origin, 
then there are four operations re- 
quired by the transform; 

1. TRANSLATE world to screen 
origin (by -eyex, -eyey, -eyez). 

2. ROTATE anti-clockwise about 
screen y-axis (T). 

3. ROTATE anti-clockwise about 
screen x-axis (90-9). 

4. CONVERT from left-hand toright- 
hand (change sign on all z). 

This process shows the need for 
mixed coordinate representations 


Fig 5 


es a aa projected point 


— the eye point is needed both in cartesian world 
coordinates (eyex, eyey, eyez) and polar world 
coordinates (r, T, a). For this 


article weshallassumetheeye yin = vector(); 


approach is required. 

In this case a perspective 

projection is appropriate where 
objects that are further from the viewer 
ve appear smaller (the further you stand 


back from your house the smaller it 
appears). This is accomplished by scal- 
ing the x and y values of each vertex 
depending on the z component. 

Fig 5 shows a side view of a scene 
viewed from an eye point. The screen in 
the middle represents the viewing plane 
(i.e. the computer screen) and it can be 
seen that objects further from the eye are 
projected smaller onto the screen. Thus if the eye is 
at (0,0,0) and a given point is at (x,y,z) which is D z- 


eye 
(0,0,0) 


/* create temp storage */ 


position input as polar coordi- vout = vector(); 


nates defined in Fig 3 and the 
cartesian version is generated 
from these values using sim- 


ple trigonometry; { 
vin[0 
pa : ; vin[1 
eyex = r sin(g) sin(T) vin[2 
eyey = cos(9) 
eyez = r sin(@) cos(T) yout 


The eye transformation is 
built up as a 4x4 matrix by 
using the matrix library func- 
tions translate, rotateyand } 
rotatex ona 4x4 identity ma- 
trix (see Fig 4) to produce Ma- 
trix eye trans. 


Projection 

Before the eye transformation can be used to trans- 
form the world data structure, the projection type 
required must be considered. The concept of a 
projection comes from the cinema — the scene is 
defined by the film and projected onto the cinema 
screen. Similarly in 3D graphics, the scene is de- 
fined by the data structure and projected onto the 
computer screen. 

If the transformation is used as created then this 
produces a parallel projection where there is no 
concept of distance from the eye (i.e. z values are 
ignored). This projection is appropriate for CAD 
systems where dimensions have to be accurate, but 
in an image rendering system a more realistic 
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] 
] 
] 


t_verts [index] .x 
t_verts [index] .y 
t_verts [index] .z 


for (index = 0; index < world.vertices; index++) 
/* do for all vertices */ 


vertices [index] .x; 
vertices [index] .y; 
vertices [index] .z; 


/* temp storage */ 


vecmatmult(vin,eye_ trans); /* do transformation */ 


vout [1]; 
vout [2]; 


freevec (vin) ; 
freevec(vout) ; 


/* free up temp space */ 


units from the screen, the x and y values are scaled 
by (D/z) before being plotted, an object at infinity is 
plotted with x=0 and z=0, and an object at z=D (i.e. 
in the plane of the screen) is plotted with un- 
changed x and y. 

The perspective 
transform can be in- 1 
troduced into the 
eye transform as a 
scaling in the fourth () 
column of the ma- 
trix by multiplying (0) 
eye trans by 
Tpers: 0 


it 
0 
0 


Fig 6 


vout[0]; /* new transformed vertex array */ 
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Clipping 

It is often necessary to produce a ‘windowed’ view 
of a scene as the screen is not of sufficient 
resolution to display a complete rendered image or 
we wish to look at a magnified portion of a scene. If 
the rendering program had to process all pixels 


when such a view was required, alarge proportion Fig8 


Eye point 


Near Clip Plane 


of the time could be spent producing pixels not 
visible. More efficient use of the renderer can be 
made if the invisible data is removed prior to 
display, a process known as clipping. Fig 9 
What is the viewing volume? 
In general we define a viewing 
volume which defines the 3D 2 
space in which we are interested; Case1 ~—~_ 
any detail outside the volume is to 
be discarded before invoking the 
renderer. The view volume is 
defined by the size of the 
viewport, or screen, the type of 
projection and user definable near 


2D WINDOW 


Case 2 


Line totally outside 


Vertex Transformation 

Once the viewing transformation has been con- 
structed it is a simple case of multiplying all world 
vertices by the transform matrix to generate a new 
set of (screen) vertices (Fig 6). The new vertices are 
held in another data structure, t_verts, which is 
a parallel of the vertices array, so it can be 
referenced by the ‘polygons’ structure. This means 
the higher levels ofrepresentation (obj ect's,poly- 
gons) do not change — only the vertices array is 
transformed. 


Clipping 

Before the transformed vertices can be displayed, 
the issue of clipping must be addressed. Given an 
arbitrary scene there are bound to be parts of it 
which the eye cannot see (they are too far left or 


discussion of the many techniques applicable, see the 
References (Part Four, October issue). 


Edge level clipping 
Consider the 2D window in Fig 9. The intersection of a 
line with the window can only be in one of the forms 
shown, any clipping algorithm must 
cope with all cases. For a simple 
rectangular window the cases can be 
identified with coordinate comparisons. 
When a line intersects with the 
window an intersection calculation 
must be performed to determine the 
line segment inside the window. This is 
computationally expensive and 
numerous imaginative techniques have 
been developed which tackle the 
problem. Several, based on recursive 
line subdivision, are suitable for 
implementation in hardware. 

\ The simple 2D case can be extended, 
Far Clip Plane greatly increasing complexity, to cope 
with a 3D parallelepiped, such as the 
view volume defined by oblique projection. The 3D 
frustrum form of the viewing volume produced by a 
perspective projection, see Fig 8, can be clipped to 
using a three stage process. Image space is distorted 
to transform 
the frustrum into a 
parallelepiped, the clipping 
is performed in the same 
manner as for oblique 
projection and the inverse 
distortion is applied. In this 
way, relatively simple 
efficient clipping algorithms 
can be made to cope with 
perspective projections. 


and far clipping planes (Fig 8). 
Several levels of clipping exist. 
Our renderer could make use of 
object, polygon or edge level 
clipping, the edge level being 
dissection of individual polygon 
edges which intersect with the 
viewing volume boundary planes. 
We shall provide a brief insight 
into the implementation of edge 
level clipping; for a detailed 


Case 4 


\ Case 3 


\ Line partially inside 


iN 


Line Pier « 


Line partially inside 
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Our rendering program is 
designed to produce a solid 
illuminated view of a simple 
scene. Once transformed to 
image space coordinates the 
scene is constant-aspect-ratio 
scaled to provide the best fit 
of the viewport, obviating the 
need for complex low level 
clipping algorithms. 
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for (object = 0; object < world.objects; object ++) 
for (polygon = objects[object] .fromp; 
polygon < objects[object].top; polygon++) 


{ 
vO = t_verts [polygons [polygon] .v0]; 
vl = t_verts [polygons [polygon] .v1]; 
v2 = t_verts [polygons [polygon] .v2]; 
v3 = t_verts [polygons [polygon] .v3]; 


plot_line(v0.x, v0.y, v1.x, vil.y); 

plot_line(v1.x, vl.y, v2.x, v2.y); 

plot_line(v2.x, v2.y, v3.x, v3.y); 

plot_line(v3.x, v3.y, v0.x, v0.y); 
/* draw polygon */ 


Fig 7 right, too far away or maybe even behind the eye). 
The process of removing these points from those to 
be plotted is called clipping (see panel, ‘Clipping’, 
page 280). 

For this article we will assume the scene is 
always completely visible and that the only clip- 
ping required will be taken care of by the graphics 
hardware (the hardware should ignore a request to 
plot a point that is off the screen). 


Wireframe 

Once all the vertices have been eye transformed 
we have a data structure which we 
can display on the 
screen. If we simply 
plotthe xand ycom- 
ponents of all the 
transformed vertices 
a scene will be dis- 
played made up of 
dots. This can be 
turned into a 
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Twice the capacity of 
your existing hard disk. 


Whatever the size of your hard disk, 
simply installing Stacker will 


Stacker. Double the capacity of 
your hard disk — in moments. 


So, installing Stacker means you can 
suddenly find room for all those 
applications you always promised yourself. 
Without the cost of buying a bigger disk. 


double its storage size. xj, Data compression 

So a 20MB disk becomes - 

without data 

40MB, 40MB becomes 80MB, : 
80MB becomes 160MB, and HM corruption. 
so on. Without the need to Using its patented LZS data 
reformat your disk. compression technology, Stacker 

In fact, PC Digest’s 5-Star pis ; removes redundancy and 

t s data is stored on disk, ; 

rating* recommended Stacker Stacker instantly compresses Saves space. And no data is 
as the easiest to learn and phlariohtradriia poland lost in the process. 
easiest to use. With the fastest expands it. You don't Just as importantly, 
performance and lowest ‘dao! eclaciies Stacker works quietly and 
memory requirements. tly in the 

It means Stacker won’t make a big background, compressing all your files. 
issue out of simply doubling your available Then, when you need them, expanding your 


_ storage space. Which is a huge benefit 
_ when you consider the memory 
athe of the latest sophisticated 


data back to normal size. Compressing it 
again, as it goes back to disk. 


STACKER: 
SOFTWARE 


DOUBLES 


In addition to the software-only version, Stacker is 
now available in two co-processor-with-software 
versions: Stacker AT/16 and Stacker MC/16, each 
featuring the unique Stacker LZS chip. 


Quickly. Without interrupting your S g 
work flow. And safely, because its in- 
built 100% compatibility stretches right 
across the board. Including the latest 
versions of DOS and Windows, of 
course. Where the need for 
compression really counts. And speed 
is of the essence. 

And if your PC hasn’t a spare slot, 
or you’re using a portable, laptop or 
notebook, the Stacker software solution — 
is just what you need. On the other — 
hand, if your aim is for a smaller driver | 


size and the highest possible speed and — 
compression ratio, you'll want to install 


-one of our two co-processor-wi 


wireframe (Fig 7) by joining all the dots (assuming 
the presence of a line draw function — 
DlLoe Linetiromx, “fromuy, bps x;-Lovy) 


Summary 

In this article we have moved from our basic data 
structure through the viewing transform to gener- 
ate a wireframe. The topics of matrix manipulation 


Matrix Transforms 

This section provides details of the form of the 
homogenous transform matrices to rotate, scale 
and translate 3D world space vertices. 


Rotation matrices 

Matrices to rotate an object about any of the 
three coordinate axes are required. They all use 
the upper left 3x3 sub-matrix and thus can be 
used in homogenous and non-homogenous 
representations. 

To rotate about an angle 7 the following 
matrices would be used for the x, yand z axes 
respectively: 


Translation matrices 

These matrices are used to provide rotation 
about an arbitrary vertex or to shift the entire 
scene to a different origin. The example matrices 
show the format for making the origin of the 
scene, T, and how to revert back to the original 
origin, respectively. 


Scaling matrices 

The leading diagonal of the transformation 
provides scaling terms for the x, yand z axes 
and a global scaling term, from top left to bottom 
right respectively. 

The example matrices show general scaling in 
each of the coordinate axes and globally, where 
the second matrix S is a reciprocal when 
compared to the terms in the first matrix. 


Matrix Library 


3D GRAPHICS PROGRAMMING PART 2 


and clipping have been introduced and applied as 
required. 

Next month we start to add the special effects to 
create a realistic image. 


A complete Reference Guide to accompany this 
series will appear in the fourth and final instalment 
in the October issue. 
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Scaling 


It is convenient to define a library of 3D matrix manipulation routines as these are frequently used in 
3D graphics. These are not described here explicitly, but the functions required are: 


Vector vector(); /* reserve space for a 4x1 vector */ 
Matrix matrix(); /* reserve space for a 4x4 matrix */ 
void freevec(); /* free space */ 
void freemat(); /* free space */ 


Matrix matmatmult (); /* multiply two 4x4 matrices */ 
Vector vecmatmult (); /* multiply 4x1 by 4x4 */ 


Matrix rotatex(); /* rotation about x-axis */ 
Matrix rotatey(); /* rotation about y-axis */ 
Matrix rotatez(); /* rotation about z-axis */ 
Matrix translate(); /* translation */ 


These assume the use of two new types: 


typedef float * Vector; 
typedef float ** Matrix; 


Note the use of homogeneous coordinates. A point in 3D space is represented by a 4x1 vector and 


a 3D transformation is represented by a 4x4 matrix. 
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